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PSL Icing Session 


Turbofan Ice Crystal Rollback Investigation and Preparations 
2:00 PM_ Leading to the Second, Heavily Instrumented, Ice Crystal Engine R. Goodwin 
Test at NASA PSL-3 Test Facility 


Determination of engine recovery after an ice accretion rollback, 
2:30 PM _ engine performance deterioration and health monitoring using D. Walker 
minimal instrumentation during icing testing at NASA Glenn PSL-3 


Preliminary Results from a Heavily Instrumented Engine Ice Crystal 


Icing Test in a Ground Based Altitude Test Facility A. Flegel 


Modeling of Highly Instrumented Honeywell Turbofan Engine 
Tested with Ice Crystal Ingestion in the NASA PSL 


Test Point Selection for Engine Crystal Icing Test at NASA PSL-3 
for focused sensitivity, peak intensity, and anti-ice evaluations 


J. Veres 


D. Dischinger 


NASA Glenn Propulsion Systems Lab Ice Crystal Cloud 


Characterization Update (2015) J. Van Zante 


Aircraft engine icing instrumentation used in the NASA Glenn 
Propulsion Systems Laboratory 


5:00 PM 


T. Bencic 
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Outline 


- Facility Description 

- Cloud Characterization Parameter Space 
- Aspects Characterized 

- Conclusions 
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PSL Facility 
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Configuration 1: 


* Engine 
* 27:1CR 


TRANSITION 


Configuration 2: 
STEAM + Component 
(FULLY MIXED) . 1110R 


* 22:1 CR with bulletnose* 
* approx, representation 
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PSL Parameter Space 


Airflow Conditions Spray Conditions 
- (Duct Geometry) - Nozzle Type & #: Mod1, Std 
- Pressure Altitude, PO - Water Pressure, Pwat 


¢« Mach, Air Mass Flow Rate, Wa - Air Pressure, Pair 


- Temperature, TPL - Air/Water Temp, Tair, Twat 

- Relative Humidity, RHPL - Water Source: City, DI 

PSL is Isentropic & Adiabatic - Spraybar Cooling Air and 
Pressure 


Physics of the Process: 

¢ Liquid water issues from the spraybars. 

¢ Water particles immediately start to evaporate. 

¢ Particles start to chill/freeze as they travel through the plenum and into the 
contraction. 

¢ The vapor... 
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Setting Conditions 


Customer Provides NASA Actions 

Ranges at Plane of Interest - Calculate Wa to set Facility 

(inlet to Fan or Booster): Conditions 

¢ Ps, static pressure - Characterize Spray Conditions 

¢ Ts, static temperature - Phase of particles at Sta 1 

« Wack primarily a function of wet bulb 
. Temp, Two. New model to 

° Altitude predict (T. Bartkus). 


¢ TWC, Total Water Content 
¢ MVD (Particle Size) 
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PSL Cloud Characterization Envelops 


National Aeronautics an 


Ambient Temperature (C) 


Altitude (kft) Two regions of 
0 10 20 30 40 50 conditions simulated 
0 ———_—=_ 1 
; | app - Engine Fan Face 
— =PSL-3 Icing Op Env. (Honeywell Engine 
10 | © Cal Regions Tests LFO1 & LETT) 
| | 
° | - Aft of Fan / Inside 
al | Booster 
Oo | (Fundamental Ice 
50; “eee ----- -- Crystal Icing Study, 
-60 - Struk et. al) 
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PSL Cloud Characterization Process 


- Cloud Uniformity 
- Total Water Content 

¢ Measurements in Center 

¢ Bulk average in Cross-Section 
- Particle Size 


- Particle Phase and Temperature - T. Bencic will 
cover next talk. 


Cloud Uniformity 


AIAA’s 

2014 Van Zante & 
2015 Bencic et.al: 
Good agreement 
between these 
three methods. 


Tomography is 
the diagnostic of 
choice. 


Grid Laser Sheet Tomography 
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Tomography — near real-time monitoring 


E34 CF=2.06 E37 CF=1.78 


Laser li ght 


1.10 1.20 


Procedure: ane 
* Measure light extinction with cloud OFF (baseline) . | ma 
¢ Measure light extinction with cloud ON E dees 

(extinction due to size and number of particles) : a 
* Intensity Ratio, |;, output at every ‘pixel’ (i, j) soma 


* Calculate avg Intensity Ratio over 1x1-in Center, log 
* Calculate Concentration Factor, CF, loo/) | 
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Total Water Content: _ Measurements 


SEA, Inc. 


Multi-Wire Robust Probe lso-Kinetic Probe 
(MW) (RP) (IKP2) 


All measurements at Duct Center, TWC,, 
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MW and RP Notes 


- Meas. power to evap. impinging water 
- Minimize flow angularity into MW head 


- Correct for Collection (Collision) Efficiency, MW 
only. Effects of particle bouncing, splashing not 
accounted for. 


- Track default TWC (100% liquid) 
Added iWC (100% ice crystal) 


Phase Change Energy Requirements 
K * Pwet 


i = <$<—<$<—< Ee 
[Cig _—_ = Fisaiy.) . Leoan| *VTAS *L *W 


5100 +--BailingRaint@. 
K «* Pigae £ Sea Level 


We = 
live (To = Tomb) 7 Lrus a Ciig = = Tag) a Lovws| *VTAS *L *W 


- IWC / TWC ~ 0.88 


SEA, Inc. WCM-2000 User’s Manual 
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IKP Notes 


lso-kinetically ingests air and cloud particles (no gain or 
loss of mass) 


Evaporates all cloud particles, regardless of phase; 
measures total water vapor 


Independent measure of background water vapor is 
subtracted so that only Ice + Liquid phases are 
calculated. 


- At PSL, several sources of water vapor measurements 
at Sta 1. Largest correction applied. 


- At PSL, a radial profile of water vapor was observed 
(Fundamental Ice Crystal Icing Study, P. Struk, 3/2016) 
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sample TWC time traces 


2D 


MW (g/m3) and Ratios 
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National Aeronautics and Space Administration 
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2:01:00 


MW Response with RH sweep 
31.6 
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Total Water Content: Bulk Water Flow 


Simple Calculation based on injected Water Flow Rate, Wf, 
Air Mass Flow Rate, Wa, and Sta 1 statics. Constant C 
includes density of water. Assume uniform distribution 
across Sta 1 Duct. 


Wf (gal/min) = #Noz * Cfn * V(DeltaP) 
TWC_W+£ (g/m?) = C * Wf * Ps1 / (Wa * Ts1) 


lf add a cloud Displacement or Boundary Layer thickness 6 
TWC_W£_BL (g/m?) = C * Wf * Ps1 / (Wa * ((R-5)/R))2 * Ts1) 


TWC_WFf is the basis function for TWC measurements 
TWC_Wf_BL was recommended for LFO1 and LF11 Tests 
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Total Water Content: Bulk Meas 


Combine Measured TWC,, and Tomography CF. 
TWC_Bulk_Meas (g/m’) = ¥ (I; * (TWCgo/loo) * Ay) / > Ay 


- Created a CF curve fit based upon Pair and TWC_ Wf, 
TWC_ Bulk Fit 


IKP TWC_Bulk_Fit RP TWC_Bulk_Fit MW TWC_Em_Bulk_Fit 

7.0 5 7.0 7.0 

6.0 ° 6.0 ye tsi 6.0 y = 1.1437x 
_ . _ R? = 0.8729 R? = 0.8847 

5.0 a 5.0 5.0 
$ 0. E 2 
= 40 O = 40 Ans B 40 
iL ° x wh 2 
x! 30 O- x! 30 Bo BK “30 a 
=] oto ° = eae x yee bas 
7 20 o- 2 20 oe = a 20 kee 
2 - 9 Oo Y We: A I eS 
= 4 6000 y = 2.3897x = aa ax® & Ss 45 ae i 

: so ° R? = 0.8539 ; f® zu to 
a“ = ba 
0.0 0.0 = 00 
0 05 1 15 2 25 3 0 05 1 15 2 25 3 0 0.5 1 15 2 25 3 
TWC_WF (g/m3) TWC_WF (g/m3) TWC_WF (g/m3) 
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TWC & radial distribution of particle size 


- Scatter in data due to radial MVD effect. CFD predicts larger particles 
concentrated at center, while smaller particles more uniformly 
distributed. 


- This radial MVD profile is not currently incorporated 


Distribution of 
20 um 


100 um 
into the tomography intensity ratios. 
IKP TWC_Bulk_Fit IKP TWC_Bulk_Fit (MVD groups) 
7.0 5 7.0 A 
6.0 ° 6.0 xX 
2 5.0 - . @ 5.0 x 
F 40 ree} Z 40 2 
ia ° ir, A 016-20 um 
x 38 i 0 i la o | 625-39 um 
a, 20 aa mi 20 © «90-105 um 
S Pe 6000 S 9090 1° 
2) Pe: a 1.0 ar A> 200 um 
0.0 | 0.0 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 05 1.0 15 2.0 25 3.0 
TWC_WE (g/m3) TWC_WE (g/m3) 
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TWC Bulk Comparison 


IKP + Tomography = Water Flow + Boundary Layer 


Bulk TWC Bulk TWC 
TWC_ Bulk Meas TWC_Wf BL 
Bulk TWC Comparison With the same boundary layer 


thickness assumption, and 
same basis formulation, TWC_Wf 


oO IKP TWC_Bulk_Meas 


70 | x MW TWC_Bulk Meas 
2 60 (nine IKP: Closer to “true” TWC 
= F |. 
[50 
wn 
o 
= 4.0 
=~ 30 
ar ee MW: Closest link to 1998 
oS * ALF502 Flight Campaign 
ke 


00 10 20 30 40 50 60 7.0 8.0 


TWC_Wf_BL (g/m3) 
Used in LF01 and LF11 
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Particle Size: Measurements 


J 
Y 
AS / 4 
A 


t Probe Cloud Imaging Probe 


Cloud Drople 


(CDP) (CIP-GS) 
Refracted + Diffracted Light Shadowing 


All measurements at Duct Center. 
Future ability to shift off-center 
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DMT, Inc. 


CDP + CIP Particle Size Distributions 


aan Pee P || 
oe OMod1_CD 

AMod1_CIP 
1E+1 = 


Number Density (#/cm‘’3/um) 


1 10 100 


Particle Diameter (um) 


1000 


Normalized Cumulative Volume 


1 10 100 1000 


Particle Diameter (um) 


Sample PSD: For the CDP survey, Tair, Twat & TPL were ‘warm’. 


For the CIP survey, Tair, Twat & TPL were at test conditions, ‘cold’. 


Same Alt, Mach and RHPL. 
Pair = 20 psid, DeltaP = 175 psid; MVD = 80 um 
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Mod-1 MVD Curve Fit 


Mod1 MVD Curve Fit 


Mod1 MVD Goodness of Fit 


120.0 120.0 
100.0 100.0 
80.0 5 80.0 
= S 
& 
Ss => 
5 60.0 rem 
LL 
2 oO 
S <= 40.0 
40.0 O 
@® MVD(um) 
20.0 3 1:1 Line 
20.0 mY er 
0.0 
0.0 0.0 200 40.0 600 80.0 100.0 120.0 
0 50 100 150 200 250 #300 
DeltaP (psid) Measured MVD (um) 


For LF11 MVD sweeps: Given a Pair line, increase DeltaP. This 
“guarantees” an increase in MVD, even if exact values are not known. 
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Additional PSD Measurements aera cine: 


Phase Doppler Interferometer 
PDI 

¢ Particle size (liquid only) 

at ¢ Particle velocity 

f >, oo" pS a =— ¢ Number density 
——- im PO ° LWC 


- ; z 
PDI Transmitter 


High Speed Imager 

HSI 

* Particle size (ice & liquid) 
e Shape 

¢ Number density 

¢ TWC 


' 
cee CREE 


ee 


5694 
see 5594 ae =) Sa 


‘ i | , @e * @@eee- 
Both PDI and HSI are non-intrusive. 1624 gg43 6643 «8646 ° 6668 6730 673g «7536759 858g 
Have taken data in two most recent efforts @e++- e000 
at center and off-center. Will be reported in future. Sa me am om 
onen . ooo 


gsi 8801 8913 8963 
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Particle Phase and Temperature 


- Raman Spectra can evaluate bonded structure of water in 
both liquid and solid phases, as well as temperature. 


- Benchtop success 


Raman Spectra of DI Water Raman Temperature Calibration 


—w— Liquid 
24 
—— Ice 
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Particle Phase, Temperature - Raman 


- “Point” measurement at beam waist 
- Some success in PSL, with particles moving at 0.5 Mach 


- Continuing development 


Raman Water Spectra 
Indicates presence 
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Conclusions 


- Cloud Cal Space is 12-parameters with complex interactions 


- Radial variations (concentrations in center) noted in 


- TWC, total water content 
- PSD, particle size distribution 
- RH, relative humidity 


- Near-independent verification of Bulk TWC measurement: 
IKP2 + Tomography within 5% of Water Flow with Boundary 
Layer 


- Lowest confidence measurement is ice crystal PSD, MVD. 
Actively investigating alternate methods. 


- Some success with measuring Particle Phase, Temp. via 
Raman Spectroscopy 
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